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1) Introduction 

• An hyperbolic partial differential equation like the Burgers equation 

u 2 

(1) d t u + d x (F{u)) = 0, F(u) EE y 

exhibits shock waves (see e.g. |17|). es£ discontinuities propagating with finite velocity. 
In order to select the physically relevant weak solution, it is necessary to enforce the 
so-called entropy condition 

(2) d t ( V (u)) + d x (C(u)) < 

as suggested by Godunov |18| and Friedrichs and Lax [T5]. In the relation (j2J), r/(») is a 
strictly convex function and ((•) the associated entropy flux (see e.g. [17J, [10J or [25J). 
For the Burgers equation, the quadratic entropy is usually considered 

(3) V(u) = \ , C(«) = y • 

• The computation of weak solutions of hyperbolic problems with the lattice Boltzmann 
scheme as described e.g. by Lallemand and Luo [21] is a difficult task. A first tentative 
has been proposed in the pioneering work of d'Humieres The study of nonlinear 
scalar equation with the help of the lattice Boltzmann scheme has been emphasized by 
Buick at al [B] for nonlinear acoustics. The approximation of the Burgers equation with 
a quantum variant of the method has been presented by Yepez [31] and we refer to Duan 
and Liu [12J for the approximation of two-dimensional Burgers equation with the lattice 
Boltzmann scheme. The extension for gas dynamics equations is under study with e.g. 
the works of Philippi et al, [29], Nie, Shan and Chen J2E] and Karlin and Asinari |21] . 

• In this contribution, we start from the mathematical framework developed by Bouchut 
[3 J making the link between the finite volume method and kinetic models in the framework 
of the BGK [JJ approximation. The key notion is the representation of the dual entropy 
with the help of the discrete velocities of the lattice. In section 2, we recal this framework 
with emphasis to the one-dimensional scalar case and in section 3 we derive three equilibria 
for a kinetic distribution associated with the lattice Boltzmann method. In section 4, we 
recall the numerical D1Q3 scheme and in section 5 the link between this scheme and the 
finite volume approach. We present our numerical experiments with the Burgers equation 
in section 6 and propose natural extensions to other systems in section 7. 



2) Kinetic representation of the dual entropy 

• The Legendre-Fenchel-Moreau duality is a classic notion defined when we consider 
a convex function ??(•) of several variables. We can apply the duality transform that 
suggests that any convex function rj(») is parametrized by the slopes of the tangent 
planes. In other terms, 

(4) n*{<p) = sup (<p.W - v(W)). 
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The maximum is obtained (when it is not on the boundary of the domain of variation of 
the state W) by solving the equation of unknown W : 

(5) rf{W) = <p. 

A first example is simply rj(w) = e w at one space dimension. Then e w = if, r?*(y?) = 
ip log ip — ip and we recover in this way the fundamental tool to define the so-called 
"Shannon entropy" [30J. 

• We can derive the dual function : if dr/(W) = p»dW then dn*(ip) = dip»W and 
the "physical state" W is the jacobian of the dual entropy. In an analogous way, we can 
introduce (see e.g. [17| . jlOJ or j25]) in the context of hyperbolic conservation laws 

(6) 8 t W + d x (F(W)) = 

the so-called "dual entropy flux" C*(v 9 )- ^ * s defined with the help of the "physical flux" 
F(u) according to 

CM = <P*F(W) - C(W), 
with the condition (jSJ) as previously. Then 

(7) d(*(<p) =d<p.F(W) 

and the physical flux F(W) is the jacobian of the dual entropy flux. In other terms, all 
the physics associated with the conservation laws (jSJ) can be expressed in terms of the 
dual entropy rf and of the dual entropy flux Q*. The example of Burgers equation 
with the quadratic entropy and associated flux gives without difficulty 

(8) rf{ip) = ^ , CM = y- 

• Independently of the framework relative to hyperbolic conservation laws, the Boltz- 
mann equation with discrete velocities has been studied by Gatignol [16j (see also [7J). In 
this contribution, we write this model for J velocities in one space dimension : 

(9) dtfj + vjdjj = Qj(f) , 0<j<J. 

The unknown quantity fj(x, t) is the density of particles at point x and time t with a 
discrete velocity Vj. We have J = 2 for the D1Q3 lattice Boltzmann scheme (presented 
in section 4). The equation admits TV microscopic collision invariants M^- : 

J^MkjQAf) = 0, 1< k<N 

j 

and N — 1 for a scalar (e.g. Burgers) equation. The N first conserved moments : 

(10) W k = J2 M kjfj, l<k<N 

j 

satisfy a system of conservation laws : 

(11) d t W k + d x (^2M kj v j f j > j=0, l<k<N. 
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Of course, we make the hypothesis that this system admits a mathematical entropy r](W) 
with an associated entropy flux (,{W). We denote by ip the derivative of the entropy (id 
est drj = if • dW) . Then the following scalar expression : 

N 

(12) ip • Mj = ^f k M kj , < j < J, 

k=l 

is well defined. We call it the j° "particle component of the entropy variables". 

• The link between the Boltzmann models and the entropy variables has been first 
proposed by Perthame [28j. We follow here the approach developed by Bouchut [3]. We 
suppose that there exists J convex scalar functions h* such that 

(13) hj(<fi»Mj) = rf(ip) , v j h*j • Mj) = C(f) > V^. 
j j 

We introduce hj(fj) = sup y (y fj — h*(y)) the Legendre dual of the convex function 
h*(m). The function hj(m) is a real scalar convex function and we can write here the 
relation §5§ making for each j the link between fj and ip • Mj under the scalar form 

(14) h'j(fj) = <fMj, 0<j<J. 
The so-called microscopic entropy 

j 

is a convex function in the domain where the hj's are convex. When the hypothesis (1131) 
is satisfied, we can prove a discrete version of the Boltzmann H-theorem. If 

(15) X>*(/i) < o, 

j 

we have dissipation of the microscopic entropy : 

(16) W) + 4(^«iM/i)) < 

and this function is a natural Lyapunov function. The equilibrium distribution f^{W) 
is naturally defined by 

(17) fp(W) = (hj)\<p*Mj) , < j < J 

because the relation (J7]) holds. Then we recover the Karlin et al [22J minimization prop- 
erty : H(f) > H(f c *) for each / such that ^ . M kj f 3 = ^ M kj = W k with 
1 < k < N. 

• We recall that the equilibrium distribution of particles, the conserved variables and 
the macroscopic fluxes are given respectively by the relations f fT7|) . ffTO]) and 

F k (W) = M k jVjfp , 1 < k < N. 
j 

The macroscopic entropy and associated entropy fluxes satisfy 
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When the Boltzmann equation with discrete velocities satisfies the so-called BGK hypoth- 
esis | I . id est 

(is) QAf) = lif?-fj), o<j<j 

for some constant r > 0, the Boltzmann H-theorem is satisfied. We give the proof for 
completeness : we first have the following convexity inequality 

{Km-hm) Or-/,) >o> o<j<j. 

If the BGK hypothesis (ITS]) occurs, we have by summation over j, 

= J2 (fMj) (f? - /,) = <p.j2 M i {fP - fi) = 

3 3 

and due to (JHJ), the hypothesis (TT5|) is satisfied and the H-theorem is established in this 
case. 

• As a summary of this mathematical section in the case of the Burgers equation ([1]) 
and following the work of Bouchut j2], if there exists convex functions h*Aip) of the 
entropy variable <p such that 

(19) X>» = v*{<p) = y> E"i%) = = y 

then the equilibrium fj q (u) = defines a stable approximation in a sense detailed in 
Chen et al [8J and extended by Bouchut [1]. 



3) Particle decompositions for the Burgers equation 

• We propose in this contribution to construct kinetic decompositions of a scalar variable 
in order to solve the Burgers equation in cases where weak solutions can occur, id est 
when shock waves can be developed. We consider only the simple D1Q3 stencil with 
three discrete velocities —A, and A. Recall that the scalar A = 4| is a fundamental 
numerical parameter that is very often taken equal to unity by lattice Boltzmann scheme 
users (see e.g. [21] )• For the Burgers equation (CQ) a possible mathematical entropy is the 
quadratic one (JHJ) . The dual entropy r]*(<p) and the associated dual entropy flux C*(v 3 ) 
are given according to the relations (J5]). Due to the method proposed in the previous 
section, we search three convex functions h*^(ip) , h^ip) and h*_(ip) such that 

(20) h\(<p) + h*M + h*_(<p) = £, A - hl(<p)) = ^ 

• A first possible solution of the previous system consists in introducing some parameter 
a such that < a < 1. Then we consider the particular function 

2 

(21) h*M = (1 - a) ^ . 
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Due to (120|) . the two other dual fuctions h+((p) and h*_{ip) are determined : 



h*, = a — + 
+ 4 12 A 



n = a — 



(22) - + , _ : - - 4 12A 

The associated dual functions can be explicited without particular difficulty : 



(23) 



KU+) 

M/ ) 
h-(f-) 
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The three functions /i* introduced in OTi) and (122]) are convex when 
(24) | | < a A 

and the relation (124"|) can be interpreted as a Courant-Friedrichs-Lewy stability condition : 

At < — As. 

1 14 1 

The stability is in fact defined as the domain of convexity of the dual functions h* pre- 
sented algebraically by relations (12 ip ( 122]) and illustrated in Figure 1. The explicit deter- 
mination of the equilibrium distribution is then a consequence of the relation (JT7J) taking 
also into account that tp = u for the quadratic entropy. We have 
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Figure 1. Kinetic decomposition (12T|) (T22]) for the Burgers equation with a "centered" 
D1Q3 scheme (a = 5). 

• An other solution of the previous system (1201) can be obtained as follows. Derive the 
two relations in ( 120]) two times. Then 



(h* + )"(cp) = (htY( V ) + I , (K)"(<p) + 2 (^J'V) 



A 
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In order to have a better stability property than the condition ( |24|) obtained previously, 
we try to enforce the convexity condition (h*) (ip) > if | ip \ < A instead of For 
ip < 0, we propose to replace the inequality (h\)" (ip) = (/i*) (<£>) + f > by an equality 
Then (h*_)"(<p) = -f if <p < 0. We deduce (h%)"(<p) = and (h£)"(<p) = 1 + f if 
(p < 0. With analogous arguments, we obtain (h+) (ip) = (/ip) (<p) — 1 — f and 
(/il)"(<^) = when p > A. 
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Figure 2. Kinetic decomposition for Burgers equation, equilibria (l26|) for the lattice 
Boltzmann upwind scheme D1Q3. 



We construct in this way an "upwind" distribution for the decomposition of the dual 
entropy: 



(26) h* + (<p) 



6A 
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It is presented in Figure 2. The associated equilibrium distribution (1TTI) takes the form 
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By considering the Legendre duals of the relations ( |26|) . we have 

h+(f+) = luV^f- 



(28) 



ho(fo) 

h-(f-) 



A* 
3 J 



I f | \ 3 / 2 Ifol 

A / A 



-2A/_ 



with /+ > 
with / Gl 
with /_ < . 



• We observe that if a = 1 for the "centered" equilibrium for D1Q3 Burgers scheme, 
the null velocity does not contribute to the equilibrium because ho(ip) = ; this vertex of 
null velocity is no more active. In that case, we obtain a D1Q2 centered lattice Boltzmann 
scheme for Burgers equation. Then 



KM = t 



h* = K 



^ '" +VTV 4 ' 12 A' "- I 12 A 

These two functions represented in Figure 3 are convex if 

(30) M < A 

and the associated Courant-Friedrichs-Lewy stability condition states as follows 

1 * 
At < — Ax. 

I u I 

The dual equilibrium entropy function defined at relations (12T?|) are represented on Fig- 
ure 3. The associated components h + (f + ) and /&_(/_) of the microscopic entropy follow 
from (123]) in the particular case a = 1. Observe that ho(fo) is no more defined which is 
coherent with a choice of a "D1Q2" lattice Boltzmann scheme. The associated equilibrium 
particle distribution is obtained according to 

,2 
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Figure 3. Kinetic decomposition for Burgers ; D1Q2 centered scheme. 
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4) D1Q3 lattice Boltzmann scheme 

• As developed in the preceding section, we here consider three examples of stable 
equilibria in the context of the lattice Boltzmann scheme. More precisely, following the 
approach proposed by d'Humieres [11] , we discretize in space and time the Boltzmann 
equation with discrete velocities Q in the following way We introduce a matrix M that 
links particle densities fj (j = — , 0, +) and moments m^. For the simple D1Q3 lattice 
Boltzmann scheme, we obtain 

/I 1 1\ 

(32) m = M.f, M= -A A , u = f_ x + f + h = m x . 

\A 2 A 2 / 

• The first equilibrium fl25l) can be translated in terms of moments under the form 

m eq ' x = [u , — , a A 2 uj . 

When using the "upwind" equilibrium (j27p . we obtain an other possible value for moments 
at equilibrium : , , + 

m eq ' 2 = [u , — , A sgn(w) -) . 

The simpler scheme D1Q2 corresponds to the first equilibrium (125]) with the particular 
value a = 1 as proposed in relations f )3ip . We have only two components in this case : 

• The relaxation step is nonlinear and local in space : 

(33) ml = m c ^ = u , m* k = m^ + [m]? — rrik) for k > 2, 

with S2 = S3 = 1.7 in our simulations unless otherwise stated. The particle distribution 
/* after relaxation is obtained by inversion of relation f l32|) : /* = M _1 » m*. 

• The time iteration of the scheme follows the characteristic directions of velocity Vj : 

f s (x,t + At) = f;(x-v j At, t). 
This advection step is linear and associates the node x with its neighbours. 



5) Comparison with finite volumes 

• In [13] we have observed that a one-dimensional lattice Boltzmann scheme can be 
interpreted with the help of finite volumes. In the case considered here, we have 
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with a numerical flux ip(x + 4p, t\ at the interface between the vertices x and 
defined according to 
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We observe that the resulting lattice Boltzmann scheme is not a traditional finite volume 
scheme (in the sense proposed e.g. in [10j) if (S2, S3) 7^ (1, 1) because the distribution 
of particles after collision /* is also a function of the two (or one in the D1Q2 scheme) 
other nonconserved moments m 2 and m 3 as described in relations ( I33p . Nevertheless, 
if S2 — S3 — 1, we can give an interpretation of the associated flux (I34|) because in this 
case, f* = fj q for all j. 

• We observe that we can also decompose the "physical" flux F(») (see the relation (JT]) 
or in all generality) under the form F(u) = F + (u) + F_(u) with 

(35) F + (u) = \f?(u), F_(«) = -A /?(«). 

We have F+(u(x, t)) + F^(u(x + Ax, t)) = A (/+ q (n(x, t)) - fl q (u(x + Ax, t))) and 
when S2 = S3 = 1 the numerical flux ip introduced in f[3^|) admits the classical so-called 
flux splitting form : 

(36) ^p[x + tj = F+(u(x, t)) +F_(u(x + Ax, t)) . 

With this above link between fluxes and particle distributions (136 j) it is natural to re- 
interpret with classical flux decompositions as (l35p those proposed in this contribution at 
relations (125|) . fl27|) and fl3T|) . As remarked by Bouchut [5], the relations (125|) and fl3T|) are 
associated with two variants of the Lax-Friedrichs scheme (see e.g. Lax [25J) whereas the 
upwind scheme (127|) corresponds exactly to the Engquist-Osher [14J scheme ! 




Figure 4. 

(JSZD at t = 
velocity cr = 



A converging schock wave for the Burgers equation. The decreasing profile 
leads to an admissible discontinuity at t = 1. Then a shock wave with 
I develops. 
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6) Test cases 

We test the previous numerical schemes for two classical problems : a converging shock 
wave and the Riemann problem. We use the three variants (125|) . (127]) and ( 13T|) of the 
lattice Boltzmann scheme for each physical problem. 

• The first test case concerns a converging shock wave and is displayed in Figure 4. At 
time t = the initial profile Uq(x) is given according to 

( 1 if x < 
(37) u (x) = < l-x if < x < 1 

{ if i>1. 

When t < 1 the profile u(x, t) remains a continuous function of space x but when t > 1 
a shock wave with velocity o = | is present (see e.g. [17], [10] or [25]). It is a challenge 
if a lattice Boltzmann scheme is able to capture in a systematic way such a discontinuous 
solution. 
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Figure 5. Burgers equation. Instable D1Q3 lattice Boltzmann simulation for a con- 
verging shock with equilibrium (125]) associated to the parameters a = |, S2 = S3 = 1.7 
and A = 1.8. Computed values are displayed every 10 time steps. 

• The first experiment (see Figure 5) concerns the first centered scheme (125]) and the 
choice a = | and A = 1.8 for the numerical parameters. The result is catastrophic, as 
depicted on Figure 5. The scheme is unstable and diverges within a very little time after 
the solution becomes discontinuous. The reason is simple a posteriori. Observe that for 
the previous test case a = | and particular values u(x, t) > 1 have to be considered. 
But the convexity-stability condition (124"]) reads as \u\< -| and is incompatible with the 
chosen numerical values because we take A = 1.8 in the numerical simulation. We observe 
that under conditions that violate the inequality (1241) . the lattice Boltzmann scheme 
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is unstable in this strongly nonlinear situation, even if we respect the linear stability 
condition 

(38) < Sj < 2 

proposed initially by Henon |2U| . 
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Figure 6. Burgers equation. Stable D1Q3 lattice Boltzmann simulation for a con- 
verging shock with equilibrium (l25|) associated to the parameters a — ~, A = 3 and 
^2 = s 3 = 1-7- Computed values are displayed every 10 time steps. 




0.5 1 1.5 2 2.5 3 



Figure 7. Burgers equation. Same experiment as in Figure 6 except that s 2 — S3 = 2. 
The numerical algorithm remains stable but is no longer convergent. Computed values 
are displayed every 10 time steps. 
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• We repeat the same numerical experiment with a smaller time step. We take A = 3 
in a second experiment. The condition (12"4"1) is now satisfied and the scheme is stable. The 
results are correct and are presented in Figure 6. The shock is spread on 4 to 5 mesh 
points and we observe simply an overshoot at the location of the shock wave. With the 
extreme set of values S2 = S3 = 2 (if we refer to relation (138]) ). the numerical experiment 
does not give correct results because no entropy is dissipated. But the scheme remains 
stable as presented on Figure 7 : the numerical values remain inside an interval [—0.4, 1.7] 
relatively close to the set [0, 1] of correct values for this particular problem. The nonlinear 
stability condition enters into competition with the linear stability condition flHSl) . 

• With the same initial condition (137|) . we use the D1Q3 upwind version (j27|) of the 
lattice Boltzmann scheme. Now the stability condition is not so severe as in the previous 
case and we take A = 1.1. The results, presented in Figure 8, are qualitatively analogous 
to the previous one (see Figure 6). We observe on Figure 8 an alternance of monotonic 
and over or undershooting discrete shock profiles. 



1.2 - 




2.1 2.2 2.3 2.4 2.5 2.6 

space 



Figure 8. Burgers equation. Stable D1Q3 lattice Boltzmann simulation for a converg- 
ing shock with upwind equilibrium 02 7p with A = 1.1 and s 2 — S3 = 1.7. Eight consecutive 
discrete time steps. 
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Figure 9. Burgers equation. Stable D1Q2 lattice Boltzmann simulation for a converg- 
ing shock with quilibrium (p?Tj) . A = 1.5 and s 2 = 1.7. Computed values are displayed 
every 10 time steps. 



• With the same decreasing initial condition fl37|) . using the D1Q2 version (|3Tj) leads 
to results presented on Figure 9. We observe only an over-shooting at the discrete shock 
profile without any under-overshooting. 

• In a second set of experiments, we use the very simple "two steps" or "Riemann" initial 
condition. The first one is simply 

1 if x < 0.2 
if x > 0.2. 



(39) uo(ar) 



The entropic solution of this Riemann problem composed by the Burgers equation ([T]) 
associated with the initial condition (139|) is a discontinuity propagating at the velocity 
a = | (see e.g. [IT], |10] or |25|). With the numerical schemes introduced previously, 
this entropy satisfying solution is captured with a precision comparable to finite-volume 
type methods. The results are presented on Figure 10 for numerical schemes (125 p . ( 12 7p 
and (|3Tj) . On Figure 11, a zoom of the previous data shows that this moving shock is 
captured by a stencil of four to five mesh points. 
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Figure 10. The Riemann problem for the Burgers equation associated with the initial 
condition (139]) develops a shock wave. The figures shows the numerical solutions with 
the three variants of the scheme after 100 discrete time steps and parameters A = 3 and 
s 2 = s 3 = 1.7. 
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Figure 11. Zoom of Figure 10 around the location of the shock wave. 
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Figure 12. The Riemann problem for the Burgers equation associated with the initial 
condition ( 140 p develops a rarefaction wave. Numerical solutions with the three variants 
of the lattice Boltzmann scheme after 100 discrete time steps and parameters A = 3, 
s 2 = s 3 = 1.7. 

• We reverse the values and 1 in the initial condition (139|) and obtain in this way a 
new initial condition : 



(40) 



u (x) 



if 
if 



x < 0.2 
x > 0.2 



The entropic solution of (TTj) (HI)j) is a rarefaction wave : a continuous solution with two 
constant states and a self-similar component as detailed e.g. [T7], [ID] or Without 
any modification of the scheme, the numerical solution with the three previous variants 
are presented on Figure 12. At the tricky zones of the foot (Figure 13) and the top 
(Figure 14) of the rarefaction, the slope is discontinuous and the solution of the problem 
([!]) (j40p is just continuous. We observe that the "D1Q2" version of the lattice Boltzmann 
scheme exhibits a two point discrete structure ; in some sense the little number of mesh 
points of this version (l3~lj) induces some rigidity in the discrete approximation. 

• In this section relative to test cases for unstationary solutions of the Burgers equa- 
tion, we have observed two facts. First, a convex- satisfying lattice Boltzmann scheme 
associated with a particle decomposition (12"0~|) of the dual entropy is naturally stable even 
in circumstance where the classic linear analysis is a priori in defect. A precise analysis 
of the competition between nonlinear equilibrium and over-relaxation step (133]) could be 
welcomed. Second, under the convexity condition of the h* functions of the particle de- 
composition ( I20p . we observe that the entropy condition is automatically enforced. No 
so-called rarefaction shock has never been observed with the initial condition (l4"Uj) . 
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Zoom of Figure 12 at the foot of the rarefaction. 
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Zoom of Figure 12 at the top of the rarefaction. 
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7) Systems of conservation laws 

• The extension of the previous ideas from scalar equation to hyperbolic systems is a 
difficult task. We first consider the example of one-dimensional linear acoustics to fix the 
ideas. We recall that we can write this physical model as an hyperbolic system of first 
order : 

(4!) a 'fi) +a *CsV =0 ' 
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Then a mathematical entropy is simply a quadratic form that corresponds to the physical 
energy : 

(42) V (W) = £ + ^ ■ 

The entropy variables are the gradients of the entropy (142|) relative to the conserved 
variables (p, q) and we have 

(43) <P=(pA 

\ c 

The associated entropy flux (,{W) is easy to determine and ((W) = pq. The dual 

entropy 77* (93) = <p»W — r)(W) and the dual entropy flux = <p»F(W) — ((W) can 
be evaluated without difficulty and we obtain 

(44) V*{v)=ri(W), C» = C(«0; 
all is quadratic in this system ! 

• We approach the system (HTj) with a D1Q3 lattice Boltzmann scheme. We use the 
moments m associated with the same matrix M used for the Burgers equation (see (132]) ). 
The associated particle components of the entropy variables ipmMj introduced in (11 2 p are 
given according to 

(45) tp . M + = p + , ip . M = p , (p»M- = p--£. 

c o c o 



(46) 



The identities (TTBl take now the form 

^;(^.M+)+/i*(^.M )+/il(^.M_) = r/*(<p) 
A/i;(<^.M + ) - Xh*_ (p.M_) =C*(V)- 
We search a possible solution of system ( I46p with simple quadratic functions : h^y) = 
ay 2 and h* + (y) = h*_(y) = by 2 . After some lines of algebra, the previous representation 
and the above conditions (I46|) leads to 

.* / . X Q\ c 2 ( Xq 



(47) 



The functions proposed in ()47j) are convex under the stability condition : 

I c I < A . 

This inequality means that the numerical waves go faster than the physical ones, a familiar 
interpretation of the Courant-Friedrichs-Lewy condition (see e.g. [25J). A microscopic 
entropy H(f) = h+(f + ) + /io(/b) + h-(f-) can be easily derived from (HTj) with the 
following contributors : 

K(f+) = £fl, M/o) = - 1 fl /.-(/-) = %P- 

2(1-% 
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The particle distribution /j* 1 at equilibrium is a direct consequence of relations (ITT)) and 
(1471) and we have 



In terms of moments, the relations (HH|) reduce to = c\p as used classically since 
the work |24| of Lallemand and Luo ! 



(49) 



• In the case of shallow water equations at one space dimension we can apply the 
program presented above for the linear acoustic model and try to represent the dual 
entropy with the help of a D1Q3 particle distribution. More precisely, we consider the 
one-dimensional system of conservation laws due to A. Barre de Saint Venant : 

dtp + d x q =0 

2 

d t q + d x (^- + kp^ =0 

where k > and 7 > 1 are given positive constants. We introduce velocity u, pressure p 
and sound velocity c > according to the relations 

_ u 7 2 _ IP u 7 -l 

u = - , p = k p, c = — = 7 k p' 
P P 

Then the entropy 77 and the entropy flux ( satisfy 

1 J) 

(50) r)=-pu 2 ^ -, ( = nu + pu; 

Z 7 — 1 

the entropy variables <p = [d p n, d q rj) = (a, (3) can be evaluated without difficulty : 

c 2 u 2 

a = - — , p = u . 

7 — 1 2 

The dual entropy rj* and the dual entropy flux C* can be expressed as functions of the 
entropy variables : 

With the matrix M introduced at relation (152"]) . we denote by ip + , ipo and the particle 
components of the entropy variables tp»Mj and we have 

(p + — a + X j3 , Lp = a , = « — X/3. 

The unknown convex functions h* satisfy the identities (j4*6|) and take now the form 



7 



h* + (<p + ) + h*M + = K (a + ?P) 

Xh* + ( v+ )-Xh*_(^) = K L + ^y~ p. 



(52) 



• We prove in the following that the system of equations (j52p where the unknowns are 
the convex functions h* has no solution. In order to establish this property, we introduce 
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the equilibrium distributions /J q according to (ITT]) . We differentiate the relations ( 152]) 
relatively to a and /3. We obtain 

( f?(a + \(3) + /*» + f_\a -\0) = p 

(53) I \fl\a + \p)- \fl\a- \p) = pu 

{ A 2 fl q (a + A P) + A 2 /I q (« - A p) =pu 2 +p. 

and we are supposed to determine an increasing function of only one real variable 
such that 

(54) s ,__».(,„. + ,). 

Due to the elementary calculus ^p- = 7^(7 — l)p 7 ~ 2 = (7 — 1) — , we differentiate the 
relation ( 154]) relative to p and independently relatively to u. We obtain 

(55) - p (/o q )V) + ^(^ + c 2 ) = l, -u(/ cq )'(a) + ^ = 0. 

We extract the derivative (/^'(a) from the second equation of ( 155]) and report the result 

in the first equation. We deduce 000 

w + 3 c = A 

and this relation can be correct only for exceptional values of velocity and sound velocity ! 
This impossibility is mathematically natural : it is in general not possible to represent 
a function of two variables (the right hand side of relation (154]) ) by a simple function of 
only one variable. 

• As a summary of this section, the generalization of what have been done in this 
contribution for the Burgers equation is absolutely nontrivial and mathematically impos- 
sible for the familiar nonlinear system of Saint-Venant equations. One idea is to keep the 
approach as a possible approximation of systems of conservation laws. 

8) Conclusion and perspectives 

• We first propose a summary of the algebraic work that has to be done in order 
to determine in which domain a given lattice Boltzmann scheme is stable in the sense 
proposed by Bouchut [3] . If very interesting results are computed with a very good lattice 
Boltzmann scheme in the framework proposed by d'Humieres [llj, we first determine 
the conserved variables Wk = Ylj Mkj fj ■ Then the convective fluxes follow the relation 
F a k{W) = MkjVj 1 fp. First it is necessary to have a kinetic decomposition of the 
entropy and the associated entropy flux of the type rj(W) = J2jhj(fp) and Ca{W) = 
Ylj v< j hjif'P)- Second determine the entropy variables tp = VwV{W) and the one to 
one mapping between W and ip. Third evaluate the Legendre-Fenchel-Moreau duals 
h*(y) = supj(y/ — hj(f)) of the scalar functions hj(»). Fourth determine in which 
domain all the functions tp 1 — > h*(ip»Mj} are convex . Fifth report this domain in the 
/ space... 

• Second, we recall that in this contribution, we have applied the above procedure to 
the Burgers equation, a fundamental nonlinear scalar equation. Then nonlinear stability 
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does not reduce to a simple criterion on the "s" parameters of the lattice Boltzmann 
scheme. It remains open for as to understand why the discrete results with the lattice 
Boltzmann scheme are so well interpreted in terms of Bouchut's theory. Moreover, it is a 
natural question to know why the entropy condition is naturally enforced in the context 
of nonlinearly stable lattice Boltzmann schemes. 

• Third we observe that the situation for systems is very tricky ! Progress could result 
from t he use of a vectorial particle distribution as initially proposed by Khobalatte and 
Perthame in [23] and developed by Bouchut [2J for the kinetic finite volume approach. 
Observe that this idea has been also recognized as very useful in the lattice Boltzmann 
community for the approximation of thermal fluids, as suggested by He, Chen and Doolen 
[19J, Dellar [9] and used by Peng, Shu and Chew [27] among others. 
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